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ABSTRACT 

We have fit synthetic visibilities from 3-D (CO^BOLD + PHOENIX) and 1-D (PHOENIX, ATLAS 12) 
model stellar atmospheres of Procyon (F5 IV) to high-precision interferometric data from the VLTI 
Interferometer (K-band) and from the Mark III interferometer (500 nm and 800 nm). These data 
sets provide a test of theoretical wavelength dependent limb-darkening predictions. The work of 
AUende Prieto et al. has shown that the temperature structure from a spatially and temporally 
averaged 3-D hydrodynamical model produces significantly less limb darkening at 500 nm relative 
to the temperature structure of a 1-D MARCS model atmosphere with a standard mixing-length 
approximation for convection. Our direct fits to the interferometric data confirm this prediction. A 
1-D ATLAS 12 model with "approximate overshooting" provides the required temperature gradient. 
We show, however, that 1-D models cannot reproduce the ultraviolet spectrophotometry below 160 
nm with effective temperatures in the range constrained by the measured bolometric flux and angular 
diameter. We find that a good match to the full spectral energy distribution can be obtained with 
a composite model consisting of a weighted average of twelve 1-D model atmospheres based on the 
surface intensity distribution of a 3-D granulation simulation. We emphasize that 1-D models with 
overshooting may realistically represent the mean temperature structure of F-type stars like Procyon, 
but the same models will predict redder colors than observed because they lack the multicomponent 
temperature distribution expected for the surfaces of these stars. 

Subject headings: convection — methods: numerical — stars: atmospheres — stars: fundamental 
parameters (colors, temperatures) — stars:individual (Procyon) — techniques: in- 
terferometric 



1. INTRODUCTION 

The connection between the transport of energy in stel- 
lar atmospheres and the limb darkening of stellar pho- 
tospheres has been under investigation for nearly 100 
years in the case of the Sun. Since 1906 l|Schwarzschildl 
Il966) . the darkening of the solar limb has been inves- 
tigated to understand the roles of convection and radi- 
ation in the transport of energy through the Sun's at- 
mosphere. Assuming that the mass absorption coeffi- 
cient of the solar atmosphere was both wavelength and 
depth independent and that the angle-dependent inten- 
sity could be replaced by the mean intensity at each 
depth, Schwarzschild derived a center-to-limb intensity 
variation based on a purely radiative equilibrium tem- 
perature structure 
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where r is the optical depth, is the angle between 
the line-of-sight and the emergent intensity, and (3 is 
the linear limb-darkening (LD) coefficient. Schwarzschild 
showed this LD law to be more consistent with contem- 
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porary observations than a LD law base d on a n adia- 
batic equilibrium temperature structure. iMilnel l|1921f) 
improved on Schwarzschild's mean intensity approxima- 
tion and found that a radiative equilibrium temperature 
structure with improved flux conservation yielded a LD 
coefficient of /3 = 3/5 in better agreement with observa- 
tions. 

The studies of Schwarzschild and Milne were completed 
before hydrogen was recognized as the pr ir icipal con- 
stituent of the Sun's atmosph ere bv iPavn^ l)1925|) . be- 
fore the work of lUnsoldl 1(193(1) concerning the effects of 
hydrogen ionization on the stability of radiative equilib- 
rium against convection, and before the importance of 
the bound-free and free-free op acity o f the n eg ative ion 
of hyd rogen was recognized bv iWildll l|1939D . iPlasketd 
l) 19361) inverted solar LD observations to find a temper- 
ature structure inconsistent with radiative equilibrium, 
suggesting conve ction has a significant effect on limb 
darkening. Ikeenanl l(l938t) showed that this conclusion 
was acceptable only if the adiabatic gradient began at 
op tical dep ths considerably less than f ~ 2, contrary 
to lUnsoldf s calculations. The reasoning was that con- 
vection currents should transport little energy above the 
zone of instability and that radiative equilibrium should 
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prevail at smaller optical depths significant to limb dark- 
ening. Plaske tt's conclusion was further critiqued by 
iWoollevI (^941*) who showed that it necessitates convec- 
tive velocities large enough to blur the Franuhofer lines. 
Woolley attributed the differences between the observa- 
tions and the radiative equilibrium models to the fre- 
quency dependent opacity of the solar atmosphere, not 
yet included in the models. As model atmospheres were 
improved with the incorporation of line blanket ing and 
mixing-length convection, subsequent studies llMiinch 
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^ ^ gener- 
ally confirmed iWoollevi s suggestion and concluded that 
the effects of convectio n on limb darkening w ere subtle or 
insignificant. Recentlv. iCastelli et al.l l)1997j) have shown 
that an "approximate overshooting" modification to the 
standard one-dimensional (1-D) mixing-length convec- 
tion treatment enhances convective transport at smaller 
optical depths and provides a better fit to solar LD obser- 
vations than models without overshooting. Overshooting 
refers to the depth of convective penetration into layers 
of the atmosphere stable against convection under the 
Schwarzschild criterion. 

Photospheric convection and limb darkening are again 
the focus of some of the most exciting observational 
and theoretical studies in astrop hysics. Recen tly, space 
based photometric observations ijMatthews e t al. 200^ 
have put an upper limit on pressure mode oscillations 
excited by turbulent convection in the atmosphere of the 
F-subgiant Procyon. Precise ground based photometric 
microlensing observations are providing multi- band limb- 
darke ning measurements of solar-like stars ijAbe et al.l 
120031) and recent observations of this kind are challeng- 
ing m odel LD predictions for cool giants (Fields et al. 
|2003). Increasingly precise long baseline interferometric 
measurements at optical and near infrared wavelengths 
HWit^tkowski et a] .1120041 IC^shi et abllainl IPerriu et a,].l 
120041) are providing new tests of LD predictions. Tests 
of these predictions are important because they aid not 
only in the interpretation of interferometric data, but 
in the studies of eclipsing binaries, transiting extrasolar 
planets and their host stars. 

Along with these recent observational developments, 
3-D, tin ie dependent radi ation hydrodvnamical simu- 
lations ijStein fc Nordlu nd '199^, 'Robin son et all |2003() 
have been produced to study solar granulation. With the 
aid of such models , the study of spectra l line formation in 
solar granulation ijAsphmd et al.ll2fl04D has recently sig- 
nificantly revised the solar elemental oxygen abundance. 
These simulations model in great detail the transition 
from convection dominated transport to radiation domi- 
nated transport near the surface of the solar atmosphere. 
Self-consistent, and not employing the usual 1-D mixing- 
length theory, these simulations predict a spatially and 
temporally averaged temperature profile for the solar at- 
mosphere which samples upfiows and downdrafts across 
the model solar surface. Similar models have also 
been develope d for cooler and mo re evolved star s (e.g., 
Frevtag et all [2002: Ludwig et all 12002: Robinso neTall 
20041) . For the Sun, multi-wavelength LD observations 



are used to test the intensities predicted by these mod- 
els in addition to spectroscopic tests of the flux predic- 
tions. Qualitatively, multiwave length solar LD predic- 
tions from 3-D models (..Asplund et alJlT99^ fig. 3) and 



1-D models with overshooting ijCastelli et aLlllW^ fig. 
4) are in reasonable agreement with the observations of 
iNeckel fc Labd ('1994) . Both these studies find that stan- 
dard 1-D mixing-length theory predicts limb darkening 
which is too strong, most noticeably in the blue (500 
nm). A similar relationship between 3-D and 1-D limb- 
darkening predictions has been established for Procyon 
(jAllende Prieto et alJ200l hereafter, AP02). These pre- 
dictions for Procyon have until now not been directly 
tested by long baseline interferometric observations and 
this is the principal motivation for this paper. 

Procyon (a CMi, HR 2943, HP 61 421), an F-type 
subgiant (F5 IV-V, iGrav et all I2001D with a white 
dwarf companion, has been a prime target for both 
observational and theoretical asterqseisni ological studies 
l|Martic et alJll999HChabover et alJll999D because its vi- 
sual orbital solution and measured angular diameter pro- 
vide well-constrained estimates of its fundamental stel- 
lar parameters: effective temperature, gr avity, radius, 
and mass; most recently improved upon bv lGirard et alJ 
(|2000() and lKervella et all (|2004b(l . The model of AP02 
reproduces spectral line shifts, signatures of convective 
motion, measured from very high resolution (R~200,000) 
spectroscopy of Procyon. A comparison between this 3- 
D model and a homogeneous, hydr ostatic 1-D model at- 
mosphere from the MA RCS code (|Asplund et alJll997t 
iCustafsson et al.l 119751) . with mixing- length convection 
and the same effective temperature and gravity, shows 
that at 1 fim the normalized center-to-limb predictions 
are quite comparable, while at 450 nm the 1-D intensity 
profile is up to 20% fainter than the 3-D intensity profile 
at intermediate limb angles. This LD contrast can be 
quantified with respect to interferometric observations 
by comparing the predicted LD corrections for the two 
models. The LD correction is a model dependent scale 
factor between the angular size derived from the visibil- 
ity data assuming a uniformly bright stellar disk and the 
true wavelength independent angular size corresponding 
to a physical radius (see e.g., ©avis et al. 2000). From 
AP02, the monochromatic LD corrections at 450 nm are 
1.081 (1-D model) and 1.064 (3-D model), a difference 
of 1.6%. In this paper we test this prediction by com- 
paring both 1-D and 3-D models directly with precise 
interferometric visibility measurements obtained in two 
intermediate bands at 500 nm and 800 nm, and in the 
K-band at 2.2 /im. 

It is also interesting to compare the uncertainty in 
a typical interferometric effective temperature estimate 
(~100-200 K) with the horizontal temperature variation 
expected from 3-D models. The AP02 3-D model for 
Procyon shows a R.M.S spatial temperature variation of 
~500 K (8%) at a depth defined by a Rosseland optical 
depth of unity. As a result, while the bolometric fiux 
and the angular diameter provide a temperature which 
represents the bolometric output of the photosphere, a 
1-D atmosphere with this temperature is not expected 
to accurately represent the spectral energy distribution. 
For this reason, in addition to testing the 3-D model 
limb-darkening predictions, we also compare predictions 
of a 3-D model atmosphere to archival spectrophotome- 
try and Stromgren photometry. 

The broad wavelength coverage needed for the inter- 
ferometric test requires data from more than one in- 
terferometer. We use 500 nm and 800 nm data from 
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the Mark III interferometer ijShao et al.lll988 iD and 2.2 
/im data from the Very Large Telescope Interferometer 
VLTI), spe cifically its commissioning instrument VINCI 
Kervella et al. , 2004a) . In ij^lwe describe the observa- 
tions and present the observational data and in ^ we 
describe the model atmospheres and the computation of 
the synthetic visibilities and the synthetic spectropho- 
tometry. In 21 we show the comparison of the synthetic 
visibilities computed from the model atmospheres with 
the visibility data and archival photometric data. These 
results are discussed in !|Sland briefly summarized in 

2. OBSERVATIONS AND FUNDAMENTAL PARAMETERS 
2.1. Interferometric Observations 
2.1.1. VLTI Observations 

The European Souther n Observatory's Very L arge 
Telescope Interferometer f Clindemann et aLl '2003'1 has 
been operated on top of the Cerro Paranal, in northern 
Chile since March 2001. The new observations presented 
in this paper were obtained by combining coherently the 
light coming from the two VLTI test siderostats (0.35 m 
aperture) on the EO-Gl baseline of the VLTI (66 m in 
ground length). 

We used the VINCI beam combiner ( Kervell a et al.l 
l2003aj) equipped with a regular K band filter (A = 
2.0 — 2.4 um). Whi le the observations presented by 
IKervella et all i2004bD were limited to a maximum base- 
line of 24 m, the 22 new squared visibility measurements 
cover the 42-66 m range, significantly lower in the first 
lobe of the Procyon visibility curve. Both the 24 m and 
66m baseline data are listed in Tabled 

The raw data were processed using a wavelet based 
algorithm (Kervella ct al. 2QQ4a,) to obtain the instru- 
mental squared moduli of the coherence factors. The 
instrumental transfer function was deduced from obser- 
vations of four calibrator stars (HR1799, 6 Leo, ISMon 
and 30 Gem) whose a priori angular diameters were taken 
from the Cohen ct al. (1999) list. T heir lin e ar lim b dark- 
ening coefficients were taken from IClaretl Emh tables 
for the K band. The 24 r n basel ine data previously pub- 
lished in IKervella et al.l l)2004bfl . as well as three points 
from the 66 m baseline, were calibrated using Sirius as 
the reference sta r. It s we ll-determined angular diame- 
ter ijKervella et alJl2003bl) results in a small systematic 
uncertainty on the calibrated visibilities of Procyon. Fur- 
thermore, the Sirius-calibrated visibilities obtained on 
the longer EO-Gl baseline show no detectable deviation 
from the other calibrators. The other calibrators have 
small angular diameters compared to Procyon (Ri2 mas 
vs. more than 5 mas). This translates into small system- 
atic errors on the calibrated visibility data. These errors 
were carried through the data reduction process and are 
reported in Tabled together with the associated calibra- 
tor stars. 

Apart from Sirius, the calibrators of the VLTI- VINCI 
data are al l K g iants from the list established by 
ICohen et al.l l)1999D . For this homogeneous set of stars, 
the limb darkening correction LD/UD is always in the 
1.02-1.03 range. The LD correction for Procyon in the K 
band being just below LD/UD=1.02, it is almost iden- 
tical and this correction cancels out. For Sirius, the ef- 
fective temperature is higher ('^lOOOO K), resulting in 
a smaller LD correction of 1.012. In order to take into 



account properly this 0.7 % d ifference, we have used the 
linear LD correction fromlC^arct (2000). Considering the 
small amplitude of this correction (1.2%), the introduced 
systematic uncertainty is considered negligible compared 
to the VINCI error bar on the best estimate LD angular 
diameter of Procyon, 5.404 ± 0.031 mas, ±0.6% (see fgj. 
The infrared K band is an advantage in this matter, as 
the LD corrections are small, and do not change much 
over a broad range of effective temperatures. 

2.1.2. Mark III Observations 

The Mark III interferometric data presented 
here comprise visibility measurements at 800 nm 
ijMozurkewich et al.lll9911 fig. l(i)) from which uniform 
disk angular diameter fits have already been published 
(|Mozurkewich et all 119911 |2p03) and unpublished visi- 
bility measurements at 500 nm from D. Mozurkewich 
(2004, private communication). We list these data in 
Tables and 01 so that others may fit rn odels directly 
to these visibilities. Mozurke wich et all (|^91, 200^ 
do not provide a listing of calibrator stars. They note, 
however, that the calibrator stars are all smaller than the 
program stars and we expect the calibrator visibilities 
to be largely insensitive to model-dependent systematic 
uncertainties due to limb darkening. Seriously large sys- 
tematic errors seem to be ruled about by a comparison 
of the published Mark III uniform disk fits at 450 nm 
to the earlie r Intensity Interferometer un iform disk fits 
at 443 nm l)Hanbiirv Brown et al.l Il974^ , which agree 
within la. 

2.2. Initial Fundamental Parameters 

2.2.1. Effective Temperature 

Fundamental stellar parameters for Procyon are well 
constrained by direct measures of its bolometric flux 
and angular diameter. The bolometric flux is de- 
rived from an integration of Procyon's spectral en- 
ergy distribution (SED). Absolute spectrophotometry 
in the ultraviolet come from the International Ultra- 
violet Explorer (low dispersion, large aperture spec- 
tra SWP43428 and LWR09108 from the INES archive 
ijRodriguez-Pascual et al.|ll99 9^. effectively from 170 nm 
to 306 nm), and from the Hubble Space Telescope Space 
Telescope Imaging Spectrograph (10 one-dimensional ex- 
tracted spectra in the sequence o61501020 to o615010k0 
from 220 nm to 410 nm). For integration, the ul- 
traviolet data are binned to 1 nm resolution. Abso- 
lute spectrophotometry in the visual and near-infrared 
wavelength range (322.5 n m - 1027.5 nm, 5 nm r esolu- 
tion) has been obtained bv iGlushneva et all l|1992D . For 
the SED beyond 1 micron, the broad-band p hotometry 
(JHKLMN) from IMorel fc M agnenal^ fWm), with the 
the absolut e calibra t ion of [Johnson (1965) for JKLMN 
and BesscU fe BrettI (|1988j) for H, is fit using linear in- 
terpolation in log A — log F\ space. We have adopted a 
uniform flux uncertainty of 5% for all the spectropho- 
tometric data. The most carefully calibrated photomet- 
ric standards are not lik ely to be more accurate than 
4% in the absolute sense l|Bohlin fc G illiland 2004). On 
the basis of Procyon's distance (3.5 pc), interstellar dust 
extinction is assumed to be negligible. The integrated 
flux from 170 nm to 10200 nm is 17.82 ± 0.89 x lO'^ 
W m^^. This value is in good agreement with the value 
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18.08 ± 0.76 W derived bv lCode et all iITqtI . The 

difference in the integrated fiux between STIS and lUE 
over their common wavelength range is less than 1%. The 
integrated STIS flux is ^ 6% larger than the integrated 
iGhishneva et alJ l|1992D flux between 322 nm and 410 nm. 
Our integrated flux and the limb-darkene d angular diani- 
eter (^ld = 5.45±0.05 mas) from Kcrvc Ua et all l)2004bD 
yield Tcs — 6516 ±87 K. This val ue is in good agree ment 
with Teff = 6510 ± 150 K from ICode et all l(l97l . A 
comparison of the angular diameter uncertainty (0.9%) 
and the integrated flux uncertainty (5%) shows that the 
uncertainty in Tcff is dominated by the uncertainty in the 
spectrophotometry. 

2.2.2. Radius and Surface Gravity 

The HIPPARCOS parallax (285.93±0.88 mas, 
iPerrvman et alJ 11997!) and the angular diameter (0ld 
= 5.45±0.05 mas) yield a radiu s of 2.05±0.02 Rq. The 
orbital solution of ICirard et al.l l)2000l|) provides a mass 
for Procyon A, which is 1.42±0.04 M p, adop ting the 
HIPPARCOS parallax l|Kervella et all l2004bf) . The 
mass together with the radius yields a surface gravity, 
log(f;) = 3.95±0.02 (in cgs units), which is uncertain at 
the 5% level. 

3. MODEL ATMOSPHERES 

The construction of model atmospheres for Procyon 
and the subsequent calculation of synthetic radiation 
fields for comparison with interferometry and photom- 
etry is done three way s: (1) Stand-alone 1-D PHOENIX 
(jHauschildt et al."1999^ structures, radiation fields and 
spectra ; (2) CG^BDLD (Frcvtag et al. 2002,) 3-D struc- 
tures temporally and spatially averaged to 1-D, then read 
by PHOENIX for computation of the corresponding radia- 
tion fields; (3) Stand-alone 1-D ATLAS 12 (Kurucz 1992) 
structures, radiations fields and spectra. These models 
will be discussed in more detail below. 

3.1. Input Parameters 

The ATLAS 12 models are plane parallel and require an 
effective temperature Teff and a surface gravity, log{g), 
as input parameters. The spherical PHOENIX models ad- 
ditionally require a stellar radius because in the spherical 
case the luminosity, not the fiux, is constant. For Pro- 
cyon's photosphere this distinction is negligible, however 
the boundary radii of the spherical and plane-parallel 
models differ as discussed below. CO^BOLD models a lo- 
cal, 3-D, statistically representative volume of the stellar 
atmosphere in Cartesian geometry with periodic lateral 
boundaries. As input parameters C0''^B0LD requires the 
surface gravity and the specific entropy of material en- 
tering the computational volume through the open lower 
boundary. The entropy strongly infiuences the amount 
of heat entering the volume from below, and plays a role 
analogous to the effective temperature in the 1-D models. 

3.2. PHOENIX Models 

A grid of forty- two 1-D, spherical, hydrostatic mod- 
els from the PHOENIX general-purpose stellar and plane- 
tary atmosphere code (version 13.07, for a descripti on see 
IHauschildt et all 119991: IHauschildt fc Baronlll999|) have 
been constructed with the following parameters: effec- 
tive temperatures, T^ff = 6420 K to 6620 K, 10 K steps; 



log(5) = 3.95, radius, R = 1.46 x lO" cm; depth- 
independent micro-turbulence, ^ = 2.0 km s~^(a value 
suggested for 1-D models by AP02 from their analy- 
sis of Procyon's spectrum); mixing- length to pressure 
scale ratio, a = 0.5 and 1. 25; solar elemental abundance 
ijAnders fc Grevess^ligsgf) : LTE atomic and molecular 
line blanketing, typically 10^ atomic lines and 4 x 10^ 
molecular lines dynamically selected at run time; non- 
LTE line blanketing of H I (30 levels, 435 bound-bound 
transitions), He I (19 levels, 171 bound-bound transi- 
tions), and He II (10 levels, 45 bound-bound transitions); 
and boundary conditions: outer pressure, Pgas = 10~^ 
dynes cm~^, extinction optical depth at 1.2 /im: outer 
10~^°, inner 10^. For this grid of models the atmospheric 
structure is computed at 100 radial shells (depths) and 
the radiative transfer is computed along 100 rays tan- 
gent to these shells and 14 additional rays which impact 
the innermost shell, the so-called core-intersecting rays. 
The intersection of the rays with the outermost radial 
shell describes a center-to-limb intensity profile with 1 14 
angular points. The mixing-length theory for convec- 
tion in PHOEN IX is very s i milar to the M ihalas formu- 
lation (jL^idwig et al . 2002; Mihalas 1978) with no over- 
shooting. A second grid of 24 PHOENIX models has been 
constructed as components to composite models which 
simulate the effects of granulation on the spectral en- 
ergy distribution (see i l3.6l Tabled}. Other than Teff, 
these models have the same input parameters as the first 
grid. The non-LTE H and He aspect the stand-alone 
PHOENIX models is an insignificant factor for the angular 
diameter fits and synthetic spectrophometry presented 
here, as we have confirmed by converging purely LTE 
PHOENIX models for comparison. The H^ ion is treated 
in LTE for aU the PHOENIX models. The effects of treat- 
ing many thousands of lines out of LTE in a solar-type 
star h as been recently explored by IShort fc HauschildJ 
HqHS) ■ Such effects could be important for the synthetic 
spectrophotometry, however such massive-scale non-LTE 
models for Procyon have yet to be calculated and are be- 
yond the scope of this paper. Below we specifically refer 
to PHOENIX models "A" and "B" . For both models Teff = 
6530 K a nd l og(ff) — 3.95 (consistent with the values de- 
rived in ii2.2l) . but they have different values for a, 1.25 
and 0.5, respectively. 

3.3. CO^BOLD and PHOENIX Models 

The "conservative COde for Computation of COm- 
pressible COnvection in a BOx of L Dimensions (L = 
2,3)" (CO^BOLD) was mainly dev eloped by B. Frey- 
tag and M. Steffen (for details see iFrevtag et all [20021: 
iWedemever et all l2004|) . The CO^BOLD models (series 
codes d3gt65g40n2.80-a4-3 (grey) and d3gt65g40n3.10- 
26-3 (non-grey)) for Procyon have an average Tcfi = 
6500 K, a prescribed log{g) — 4.0, and solar abundances. 
For the non-grey model, five wavelength groups are em- 
ployed to describe the wavelength dependence of the radi- 
ation field with in a multi-group radiative transfer scheme 
l|Ludwiglll99 l. The group-averaged opacities a re bas ed 
on data provided by the ATLAS 6 code ( Kurucd l 19791) . 

Integration with the 1-D PHOENIX code provides di- 
rect opacity sampling at high spectral resolution for cal- 
culation of the radiation field based temperature struc- 
tures provided by CD^BOLD. The 1-D temperature, pres- 
sure, and depth arrays (91 depth points) from CO'''BOLD 
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are read by PHOENIX. The physical depths provided by 
CO^BDLD are relative to the layer where the average Rosse- 
land optical depth is unity, but PHOENIX requires absolute 
radii for construction of the 1-D spherical radiation field. 
The value 1.46 x 10^"'^ cm, corresponding to the physi- 
cal radius of Procyon, is therefore added to the CO^BDLD 
depths. The accuracy of this value is not critical because 
the thickness of the atmosphere is ^ 0.1% of the stellar 
radius. 

We take two approaches for the computation of the 
mean radiation field of the CO^BOLD models. (1) we com- 
pute a PHOENIX radiation field from a temporally and 
spatially averaged 3-D temperature structure. (2) we 
compute 980,000 horizontal positions from the spatial 
and temporal evolution of the granular flow (see Figure 
n]for a snapshot of the flow) and sort them into 12 groups 
according to their intensity (see Figure |21 for the his- 
togram of intensities). Next, we average the vertical slabs 
on surfaces of fixed optical depth. This results in twelve 
temperature stratifications representing the dark and in- 
creasingly brighter granulation areas. PHOENIX radiation 
fields are then computed for each of the twelve struc- 
tures. These radiation fields are then co-added, weighted 
by the surface areas of the twelve intensity groups (for 
weights see Table ^J. Theoretically, we expect approach 
(2) provides a better match to the actual radiation field 
of Procyon. This is because radiative transfer in an inho- 
mogeneous medium is non-linear and one cannot inter- 
change radiative transfer and spatial averaging exactly. 
In this respect, the twelve component model should be 
an improvement over a global spatial mean. 

We find the temporally and spatially averaged tem- 
perature structures from the grey and non-grey CO^BOLD 
models are nearly identical at the depths of continuum 
formation for the interferometric wavebands considered 
here: 500 nm, 800 nm, and 2.2 fim. As a result, the cor- 
responding PHOENIX radiation flelds from approach (1) 
yield fit angular diameters that differ by less than 0.03%. 
Much greater than the grey versus non-grey effects are 
the differences using approaches (1) and (2). Using the 
weighted average of the twelve radiation fields to sim- 
ulate the mean center-to-limb profile yields an angular 
diameter 0.4% larger at 500 nm compared to the global 
spatial mean radiation field. As we show in ^ approach 
(2) is slightly more consistent with the observed visibility 
data. Below we refer to models "CI" and "C2", models 
constructed via approachs (1) and (2), respectively. 

Tests indicate that the PHOENIX radiation field pro- 
duced from the CO'^BOLD 1-D average structure produces 
an outer boundary radiative fiux lower than expected (by 
15% in fiux, 260 K in Tcff) for an effective temperature 
of 6500 K. This is most likely due to the different opacity 
setups in the two codes. CO^BOLD uses the older ATLAS 6 
opacity setup. As a result, a spectral energy distribution 
from PHOENIX based on the mean 3-D structure will be 
cooler than expected from the CO^BOLD Toft value. This 
is not a critical problem for the interferometric compar- 
isons because of the insensistivity of the limb-darkening 
to Teff in this range (5500 K to 7500 K, see The 
radiative fiux mismatch is a problem with regards to syn- 
thetic photometry and this is discussed below in 

3.4. ATLAS 12 Models 



We also use ATLAS 12 models (|Kuruczll 19921 and R. L. 

Kurucz 2004, private communication) which include 
an "approximate overshooting" (AO) prescription for 
convective fiux tran sport in a mixing-length formalism 
l|Castelh et al.ll997l sec. 2.3), for comparison with the 1- 
D and 3-D models just discussed. In the AO formulation, 
the convective fiux extends to lower optical depths, re- 
ducing the temperature gradient relative to models with- 
out overshooting. Three ATLAS 12 models have been 
constructed: "D" , no overshooting; "F" , 50% overshoot- 
ing; and "E" full overshooting. All three models have the 
following parameters: Teff = 6530 K; log{g) = 3.95; the 
microturbulence is depth dependent increasing from 0.9 
km s^^ to 3.2 km at t he inner boundary; a = 1. 25; so- 
lar elemental abundance l)Anders fc: Grevessell989j) ; LTE 
atomic and molecular line blanketing; and boundary con- 
ditions: outer tross = 10^^, inner tross ~ 10^. After the 
models, with 72 depth points, were computed by opac- 
ity sampling, the complete spectrum was computed at 
resolving power of R=500,000 for 17 angular "/i" points 
describing the center-to-limb intensity profile. 

3.5. Synthetic Visibilities and Fit Procedure 

The computation of the synthetic visibilities from the 
model radiation fields sim ulates the ba,ndwid th-smeared 
squared visibility (see e.g., Davis et "311120001) . At a pro- 
jected baseline B and mean wavenumber A^^, the syn- 
thetic squared visibility, 



(2) 



is computed at each wavelength from a Hankel transform, 

V{B,X)^ f 5(A)/(a*,A)JoUld(S/A)(1-/z2)1/21 ^d/^, 

(3) 

where I(n, A) is the model radiation field (in photons 
cm~2 s~ sr~^), ^(A) is the instrument sensitivity curve, 
/i is the cosine of the angle between the line of sight and 
the surface normal, and ^ld is the limb-darkened angular 
size. The mean wavenumber is computed from 



/„°° A-i 5(A) Fx dX 
J,^S{X)Fx dX ' 



where 



Fx = 2n 



I{fi, A) /i d^ 



(4) 



(5) 



is the fiux. For each model radiation field, non-linear 
least squares fits are performed to the visibility data at 
500 nm, 800 nm, and 2.2 /im, using the Interactive Data 
Language (IDL) routine CURVEFIT, yielding three values 
for 6'ld- 

The baselines for the VINCI visibilities are listed in 
Table ^ The reciprocal of the mean wavenumber for 
the VINCI transmission (including: K-band filter, detec- 
tor quantum efficiency, fibers, atmosphere, and MONA 
beam combiner) is Aq — 2.182±0.002 /im using an appro- 
priate synthetic spectrum for Procyon (a Toff=6530 K, 
\og{g) = 3.95 PHOENIX model). An accurate estimate for 
the mean wavenumber is important because the fit angu- 
lar diameter scales linearly with the mean wavenumber. 
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The Mark III sensitivity curves are assumed to be 
Gaussian with central wavelengths A n = 500 nm and 800 
nm, each with a FWHM of 20 nm f Mozurkewic h et al.l 
[T99l|) . The Mark III bands are '^20 times narrower 
than the VINCI band, therefore the mean wavenumber 
is much less sensitive to the shape of Procyon's spectral 
energy distribution. 

3.5.1. Accounting for Extension in the Spherical Models 

A proper comparison of the best fit 0ld values derived 
from PHOENIX, CO^BOLD +PHDENIX, and ATLAS 12 mod- 
els requires a correction to the stand-alone PHOENIX val- 
ues because of spherical extension effects. The PHOENIX 
structures extend outward to a radial shell ~ 0.4% above 
the Tioss — 1 radius while the model plane-parallel struc- 
tures have very small extensions for a fixed gravity (the 
CO^BOLD and ATLAS 12 structures have extensions of 
~ 0.1%). The radiation field I{fi,X) refers to the an- 
gular distribution of intensities emerging from the out- 
ermost model layer. As a result, the spherical models 
with greater radial extension yield larger angular sizes 
relative to the less extended plane-parallel models. This 
is true even when the angular size of the continuum- 
forming radii for the two model geometries agree because 
their outer boundary radii differ. Thus, the PHOENIX 0ld 
values have been scaled down to correspond with the 
Tross = 1 radius (in other words reduced by 0.4%) for 
comparison with the plane-parallel models. For red gi- 
ants the correc tion is much larger, for e xample 7% in the 
case of Tp Phe ijWittkowski et alJl2004j) . This correction 
reminds us that even relatively compact atmospheres are 
"fuzz y" and that the stella r radius must be carefully de- 
fined fBa schek et"al]ll99H) . The need for this correction 
also speaks to the precision of the observational data. 

3.6. Synthetic Spectrophotometry and Photometry with 
3-D and 1-D Models 

At present, the coarseness of the wavelength-dependent 
opacity in the CO'^BOLD model does not permit direct pre- 
dictions of accurate colors. Therefore, we use the spa- 
tial intensity distribution of the 3-D model as a guide 
to construct multicomponcnt models of Procyon's gran- 
ulation pattern with PHOENIX. From the grey CO^BOLD 
model a snapshot of the emergent intensity in integrated 
("white") light is shown in Figure ^ Figure [3 depicts a 
histogram of intensities, including 12 intensity intervals 
for which we construct PHOENIX models (see Table 0] for 
corresponding surface area weights). The bimodal nature 
of the histogram shows that the model photosphere can 
be divided roughly into "cool" and "hot" components, 
as also indicated in the intensity map. In fact, we ini- 
tially experimented with such a two-component model, 
and later refined the spectral synthesis to the present 
12-component model. 

For constructing the PHOENIX models we need to spec- 
ify a Toff. Assuming a linear darkening law 



1-/3(1-m) 



(6) 



we convert emergent intensities in vertical direction to 
fluxes F according to 



F, = 7r/,(1)(1-| 



With the additional assumption of a wavelength inde- 
pendent (3 we directly obtain Toff since the grey CO^BOLD 
model provides wavelength integrated intensities. The 
Toff of the inner 10 zones are taken as the average of 
the Toff at the two adjacent boundaries. For the two 
outer zones half the interval width ATcs of the adja- 
cent zone is added (subtracted) to the hottest (coolest) 
boundary. The radiation fields of the individual compo- 
nents are added, weighted by the surface area fraction 
associated with the represented intensity interval. The 
resulting Toff values and weights for the 12 zones are 
listed in Table 01 Two values for (3 have been chosen: 
0.5 and 0.6. These values are consistent with published 
linear limb-darkening coefficients for broad-band optical 
filters (Claret 2000). More importantly, these two /3 val- 
ues yield Toff distributions which, after the twelve 1-D 
model SEDs are weighted and co-added, have bolomet- 
ric fiuxes and corresponding effective temperatures which 
bracket the CO^BOLD model (Toff = 6500 K). The com- 
posite spectrum Toff is computed as 



To, 




E 



1/4 



dX 



(8) 



where a is the Stefan-Boltzmann constant, are the 
weights given in Tabled) and TJ are the flux distributions 
from the component models. 

3.6.1. Synthetic Photometry 

Synthetic Stromgren indices llStromgrenI I1966D are 
computed using the formulae of iMatsushimal f 1 969j) . In 

this formalism, a normalization constant, or zero point, 
kij is needed to transform each natural index. 



(i - j)nat = 2.5 log' 
to a photometric index. 



FxSj{X)dX 

(9) 

FxS,{X)dX ■ 



(i - j) = {i~ i)nat + h: 



(10) 



directly comparable to observations. In the literature, 
the normalization of the theoretical indices is often based 
solely on mod els of Vega (e .g., ISmallcv & Kupka 1997), 
however iLester et al.l ()1986D have shown that these zero 
points are a function of (b — y). Since Vega and Pro- 
cyon do not have the same {b — y), using Vega leads 
to systematic errors in the synthetic indices. We also 
wish for obvious reasons not to use a model of Procyon 
to obtain these zero points. Therefore, we use two in- 
dependent S£ectro£hotom£^ of Procyon 
l|Kiehlinel IT987t IGlushneva et alJ I1992D and two differ- 
ent sets of filter sensitivity curves ijMatsushimal Il969t 
Crawford & Barnes 1970) to determine the following zero 
points: 



- y)obs - " y)„at ^ 

(Cl)obs - (Cl)„at^ 



-0.503 ±0.004 ±0.007 (11) 
-0.088 ±0.003 ± 0.009 (12) 



(7) 



These zero points are given with two sets of un- 
certainties. The first set arises from the range 
of zero points obtained using the four different fil- 
ter/spectrophotometric data combinations. The second 
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set is from to the mean observational error in the ob- 
served indices (Crawford & Barnes 1970). If we use spec- 
trophotometry of Vega l,Bohlin 1996) the mean ci zero 
point is i nstead -0.165 mag, consistent with the ci shift 
shown bv lLester et"all l|1986l ^ at the {b - y) color (0.272) 
of Procyon. 

4. RESULTS 

4.1. Multiwavelength Angular Diameter Fits 

TableOlists the best fit 6'ld values (and la uncertain- 
ties) at 500 nm, 800 nm, and 2.2 /im for seven different 
atmospheric structures and a uniform disk model. Figure 
shows all the visibility data together with the synthetic 
visibihties from the CO^BDLD + PHOENIX model "C2" with 
the mean fit angular diameter of 5.403 mas. The best fit 
angular diameters for all seven models are compared in 
Figure 01 

There is a clear trend in the fit results toward shorter 
wavelengths. The seven models yield a significantly 
wide range of angular diameters at 500 nm, differing by 
100±23 fi'AS, while differing only by 17±12 /ias at 2.2 /im. 
Of the seven models, two models stand out as best repre- 
senting the visibility data at all three wavelengths with a 
single angular diameter: CO^BDLD -I- PHOENIX model "C2" 
and ATLAS 12 model "F" (with 50% overshooting). Both 
models yield angular diameters at the 500 nm and 2.2 /im 
that are the same within the formal uncertainties: 15±17 
/ias and 4±17 /ias, respectively. Unlike model "F" , model 
"C2" has no free parameters for convection. 

4.2. Model Comparisons to Spectrophotometry and 
Photometry 

Figure I3a) shows that the synthetic SED of ATLAS 12 
model "E" fails to reproduce the observed continuum be- 
low 160 nm. This is true not only for model "E" , but all 
1-D models, regardless of convection treatment, which 
are consistent with both the measured angular diameter 
and the measured bolometric flux. While the contribu- 
tion to Procyon's bolometric flux from radiation between 
130 nm and 160 nm is negligible, the absolute flux level 
in this region provides strong evidence for a multicom- 
ponent temperature surface distribution and therefore 
granulation. 

The absolute continuum fluxes between 136 nm and 
160 nm are from the Goddard High Resolution Spec- 
trograph data sets: Z2VS0105P (PI A. Boesgaard), 
Z17X020CT, Z17X020AT, Z17X0208T (PI J. Linsky). 
These spectra were originally obtaine d to study the Pro- 
cyon's atmospheric boron abundance l|Cunha et al.l2000(l 
and chromosphere l|Wood et alJ 11993) ! The signal-to- 
noise in the continua of these spectra ranges from approx- 
imately 4:1 in the bluest data set to 10:1 in the reddest 
data set. These data are far superior to any other mea- 
surements below 160 nm because the continuum drops 
by more than a factor of 100 here, too much for the lim- 
ited dynamic range of lUE. The absolute fluxes below 
160 nm shown in Figure El were estimated by computing 
the mean flux between the emission lines in each spec- 
trum incorporating the flux uncertainties provided with 
each calibrated data set. 

The observed continuum below 160 nm is consistent 
with a spatial component of Procyon's atmosphere which 
has a color temperature in excess of the interferometric 



Tcfi'. Therefore, such a component must represent a small 
fraction of the surface, for it would otherwise produce a 
bolometric flux in excess of that observed. The contin- 
uum of Procyon between 136 nm and 160 nm should be 
photospheric, in contrast to the solar spectrum where 
the bound-free opacity of silicon causes the continuum 
to form at or slightly above the t emperature minimum, 
near the base of the chromosphere ijVernazza et al. 1976). 
Procyon's overall warmer photosphere {T^s = 6540 K vs. 
Tcfi = 5770 K) will contribute more flux to this spec- 
tral region than does the Sun. Furthermore, the warmer 
components of Procyon's photosphere, which predomi- 
nately contribute to the continuum below 160 nm, have 
a silicon bound-free opacity approximately three times 
weaker than expected for so lar photospheric conditions 
ijTravis & Matsushimalll96'M fig. 6). Thus we expect 
Procyon's continuum between 136 nm and 160 nm to 
form at depths beneath its chromosphere. 

The ultraviolet SED longward of 160 nm, where the 
multicomponent temperature effects are less prominent, 
may be impacted by non-LTE treatment of iron-group el- 
ements. iShort fc Ha uschildt ( 2005) find non-LTE models 
for the Sun have substantially (up to 20%) more near-UV 
flux relative to LTE models. Whether such differences 
would exist for a similar non-LTE model of Procyon is 
not clear. 

The lower panel of FigureElshows the ATLAS 12 model 
"E" SED matches the Composite SED B model fairly 
well in the u- and y-bands, but has significantly lower 
flux in the b- and ii-bands. In the composite SED model 
the peaks of the 12 components all contribute flux to 
the V- and 6-bands, while at longer wavelengths the flux 
contrast is lower. Table and Figure show that the 
Toff' = 6530 K ATLAS 12 models all predict (& - y) val- 
ues significantly redder than observed. In particular, the 
(6 - y) value for ATLAS 12 model "F", the best fitting 
1-D model to the interferometric data, differs by -1-0.04 
mag, significantly in excess of the estimated uncertainties 
(0.012 mag, see ^SXT|) . 

The differences in the derived indices between PHOENIX 
model "A" and ATLAS 12 model "D", particularly for 
ci, result from differences in these models' temperature 
structures (see Figure EJ. While these two 1-D mod- 
els were constructed to have the same input parame- 
ters, they differ in geometry (spherical vs. plane par- 
allel) and line formation (non-LTE versus LTE treat- 
ment of hydrogen and helium and constant vs. depth- 
dependent microturbulence) assumptions. To investi- 
gate this discrepancy, we converged two LTE, plane- 
parallel PHOENIX models with constant microturbulences 
of and 2 km s'-^. These models yield colors nearly 
identical to model "A" , suggesting that the 1-D mixing- 
length convection treatments (without overshooting) dif- 
fer between PHOENIX and ATLAS 12. A difference in the 
convection treatments of these two codes is noted by 
IShort fc Haiischildti 

5. DISCUSSION 

5.1. Model Differences, Interferometric Uncertainties, 
and Future Measurements 

As predicted by AP02, we find, based on compar- 
isons with high-precision multiwavelength interferomet- 
ric data, that 1-D model temperature structures based 
on standard mixing-length theory (without overshoot- 
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ing) produce center-to-limb intensity profiles which are 
too limb darkened at blue wavelengths relative to near-IR 
wavelengths. The magnitude of this effect appears to be 
slightly smaller than predicted, 0.4 — 1.5% versus 1.6%, 
depending on which 1-D model is adopted. Furthermore, 
limb-darkening predictions from a 3-D hydrodynamical 
model appear to be consistent with these data, also in 
line with AP02's predictions. 

Considering the three ATLAS 12 models, the fit angular 
diameters differ negligibly (0.2%) at 2.2 /im, but differ by 
1.1% at 500 nm due to different convection treatments. 
Considering all seven model atmospheres, from PHOENIX, 
ATLAS 12andC05B0LD + PHOENIX, the range of fit angu- 
lar diameters is again small (0.3%) at 2.2 /im, with differ- 
ences up to 1.9% at 500 nm. The error bars on the inter- 
ferometric data sets yield formal uncertainties of 0.1% to 
0.2% in the absolute angular diameters. It should there- 
fore be possible to choose between these models in the 
absence of significant systematic errors. Any systematic 
errors in the VLTI/ VINCI a ngular diamet ers is expected 
to be less than 0.3% (,Kervella et al.l2004bj) . The system- 
atic errors in the Mark III are not as well characterized, 
however the three 1-D models without overshooting yield 
larger angular diameters at 500 nm relative to 800 nm as 
expected, while the 3-D CO^BOLD -t- PHOENIX model "C2" 
and ATLAS 12 model "F" for 50% overshooting yield the 
same angular diameters. These relationships hold out 
to 2.2 /im, where model "C2" yields the same angular 
diameter within 1.3a. This suggests that systematic er- 
rors between the VLTI/ VINCI and Mark III data sets 
are minimal. 

Table and Figure ^ show that an angular diameter 
derived from the VLTI/ VINCI visibilities is the least sen- 
sitive to which model is employed. The mean angular 
diameter value at 2.2 /im together with la range in an- 
gular diameters (0.26%) and the upper hmit on the sys- 
tematic error in the VLTI/VINCI diameter (0.3%) yields 
a best estimate for the angular diameter of Procyon: 
5.404±0.031 mas. The corresponding best estimates for 
Teff, radius, and surface gravity are listed in Table □ 

Interferometric observations with superior calibration 
are certainly warranted to check these results. Indepen- 
dent multiwavelength angular diameter measurements 
are in order, however this is not the only approach 
presently available. Limb darkening, in addition to its 
affect on these angular diameter measurements, can be 
constrained by a precise measurement of the second lobe 
of the visibility curve. While angular diameters derived 
from the 3-D and standard 1-D models differ by only 1% 
at 500 nm, the squared visibility at the peak of the second 
lobe is predicted to differ by 10% at 500 nm. At longer 
wavelengths the predicted difference is smaller: 6% at 
800 nm and 3% at 2.2 /im. Precise second lobe measure- 
ments at 500 nm and 800 nm could possibly be made 
with the Navy Prototype Optical Interferometer (NPOI) 
and the Sydney University Stellar Interferometer (SUSI). 
Distinguishing between the 3-D CO^BOLD model and 1-D 
overshooting ATLAS 12 model "F" on the basis of second 
lobe measurements would appear to be much more chal- 
lenging. The squared visibility amplitudes at the peak 
of the 2nd lobe for these two models differ by less than 
1% at 500 nm, 800 nm, and 2.2 /im. 

5.2. Procyon's Mean Temperature Gradient 



The continuum forming region of Procyon's atmo- 
sphere is expected to have a physical thickness which 
is very small (< 0.1%) in comparison with its stellar ra- 
dius. Therefore, the wavelength-dependent uniform disk 
angular diameter is unlikely to result from an extended 
atmosphere where significantly different physical depths 
are probed at different wavelengths, as in the case of cool 
supergiants. Instead this effect must reside in the distri- 
bution of temperature with optical depth. 

The best fit angular diameter values presented in Table 
Eland Figure 0]are quite insensitive to the model effec- 
tive temperature. Fits obtained from a grid of PHOENIX 
models spanning a 2000 K range yield angular diameter 
differences of only 0.007 mas, 0.006 mas 0.004 mas at 500 
nm, 800 nm, and 2.2 /xm respectively; values comparable 
to the formal fit uncertainties. The atmosphere models 
in Table [Slyield different angular diameters because they 
have different temperature gradients, as shown in Figure 

m 

The mean temperature structure of the 3-D CO'^BOLD 
model and the 1-D ATLAS 12 model "F" have shallower 
temperature gradients in the range 1 < Ti-oss ^ 5 relative 
to the 1-D models without overshooting. This is expected 
because the radiative flux should be reduced when the 
convective flux is increased. Figure [S] shows the fraction 
of radiative flux as a function of optical depth for these 
models. The mean 3-D model and the 1-D model with 
overshooting have only ^ 80% of their total flux in radi- 
ation at Tfoss = 1. At depths Tioss ^5 8 the mean 3-D and 
the 1-D overshooting radiative flux structures are quite 
similar. These best fltting models to the interferometric 
data show significant convective flirx at Tj-oss ^5 1 we 
conclude that the interferometric data provide evidence 
for convective overshooting in Procyon's atmosphere. 

The shallower temperature gradient introduced by 
overshooting reduces the relative limb darkening between 
500 nm and 2.2 /xm. This accounts for the variation in 
the best fit 6'ld values for models with different convec- 
tion treatments. An analysis of the contribution func- 
tions, S{T\)e~'^'', at the three wavelengths shows the 500 
nm continuum forms in a range of depths partially below 
the depths forming the 800 nm continuum. The 800 nm 
band lies near the peak of the bound-free H" opacity, 
the dominant continuum opacity in the optical. At 2.2 
/im the H~ opacity, here free-free beyond the 1.6 H 
opacity minimum, is less than at 800 nm and similar to 
the 500 nm opacity. However, at 2.2 /im the intensity is 
much less sensitive to the temperature of Planckian ra- 
di ation than at 500 nm. F or the Sun, it was pointed out 
bv lGingerich et all l)197lD that the lack of a steeply ris- 
ing Planck function (with temperature) at 1.6 /xm means 
that the intensities emerging from the hot convective lay- 
ers are better probed at wavelengths immediately to the 
red of the Balmer limit (365 nm) even if the optical depth 
is a factor of two larger there than at the 1.6 /im opac- 
ity minimum. In short, similar depths are probed by 
the 500 nm and 2.2 /xm bands, but the radiation field is 
much more sensitive the temperature gradient at 500 nm 
compared to 2.2 fun. 

One-dimensional models without overshooting cannot 
be made more consistent with the multiwavelength visi- 
bility data by adjusting the mixing length parameter a. 
Based on the two stand-alone PHOENIX models "A" and 
"B" (see Table El, changing a from 1.25 to 0.5 has no 
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effect on the derived angular diameters at 800 nm and 
2.2 /im. At 500 nm the two best fit angular diameters 
differ by 0.35%. Note that model "B" with a = 0.5 yields 
a slightly larger angular diameter at 500 nm relative to 
the a — 1.25 model. Thus, reducing a slightly increases 
the degree of limb darkening for a model which is al- 
ready too limb darkened compared to the observations 
presented here. 

5.3. Photometry and 1-D Overshooting 

As discussed above, we find that the 1-D ATLAS 12 
model "F" with 50% overshooting fits the visibility data 
as well as the 3-D CD^BOLD model and that the two 
models have similar mean temperature gradients. How- 
ever, models with mixing-length theory convection and 
overshooting have been subject to extensive comparisons 
with photometry, with the results fayqring models with- 
out o vershooting ijCastelh et al.lll997t iSmallev fc Kupkal 
11997(1 . In general, 1-D convection treatments (e.g. 
mixing-length variants, turbulent convection) for model 
stellar atmospheres with otherwise identical parameters 
yield different SEDs. Stromgren photometric indices are 
often employed to facilitate the comparison of these dif- 
ferent synthet ic SEDs with a lar ge number of real SEDs. 
Such studies l|Heiter et al.ll200^ and references therein) 
generally find that 1-D atmosphere models with no over- 
shooting and a = 0.5 more closely reproduce the colors 
of real atmospheres relative to other 1-D models. This is 
demonstratecl for Procyon by ISmallev fc Kupkal l|1997f) 
who show that to reproduce the observed (6 — y) index 
with an ATLAS 9 overshooting model requires a model 
with Toff = 6830± 142 K (the uncertainly based on pho- 
tometry alone), a value in excess of their derived inter- 
ferometric value, Tpff = 65 60 ± 130 K. This confirms the 
findings of iCastelh et al.l (|W97). These results appear 
to be strengthened by our more precise value for the in- 
terferometric T^s value, now 6543±87 K, based on the 
mean K-band angular diameter value (see Tabled). Now 
the best fit 1-D overshooting model to the photometry 
would appear to require a model Tcff value even more 
incon sistent w it h the interferometric Tcff. 

Hei ter et al.l (|2002(1 find that reducing the mixing- 
length parameter a in ATLAS 9 models yields synthetic 
indices more consistent with observations. Our compar- 
isons suggest that 1-D models with lower a values better 
match photometry because this better approximates the 
composite SED in the b- and w-bands. As a result, con- 
straining the 1-D a parameter with Stromgren photome- 
try most likely constrains not the temperature structure 
of the 1-D model, but instead refiects the lack of hot- 
ter temperature components in the 1-D model, the same 
components responsible for the continuum below 160 nm. 

The composite SED models reproduce photometric in- 
dices quite close to the observed values and match fairly 
well the full UV continuum flux distribution. On the 
basis of the spatially and temporally averaged 3-D ra- 
diative flux structure (Figure IS)) one might guess a much 
redder {b—y) color for this model because of the similarly 
to ATLAS 12 model "F". However, unhke the 1-D case, 
in the 3-D case the distribution of surface temperatures 
suppresses the one-to-one connection between the mean 
photospheric SED and the mean temperature structure. 
Both composite models predict (b — y) values slightly red- 
der than observed, but much more consistent than the 



1-D overshooting models. A slightly warmer CO''^BOLD 
model (Tcff = 6550 K instead of 6500 K), more consis- 
tent with the interferometric value 6543 K, is expected to 
yield a slightly warmer distribution of temperatures and 
hence a bluer (b — y). A more consistent modeling ap- 
proach, where the colors are computed directly from the 
3-D model intensities, should clearly be pursued once the 
wavelength-dependent opacities in 3-D models improve. 

6. SUMMARY 

We find that the multiwavelength visibility measure- 
ments of Procyon obtained by the Mark HI and the 
VLTI interferometers characterize the mean tempera- 
ture structure of its atmosphere while UV and opti- 
cal spectrophotometry characterize the spatial distribu- 
tion of surface temperatures expected from granulation 
on its surface. We find that predictions from a three- 
dimensional CO^BOLD hydrodynamic model atmosphere 
of Procyon match the interferometric, spectrophotomet- 
ric, and photometric observations simultaneously, pro- 
viding confidence in 3-D model predictions which until 
now had been tested predominantly with high-resolution 
spectroscopy. Based on 3-D and 1-D model comparisons, 
the interferometric data are consistent with a tempera- 
ture structure with significant convective overshooting. 
While a 1-D model atmosphere with some overshooting 
matches the interferometric data, all 1-D models fail to 
match Procyon's spectral energy distribution, most ev- 
idently below 160 nm. Stromgren photometric indices 
appear to be sensitive to the multicomponent bright- 
ness temperature distribution on Procyon's surface. This 
would seem to severely complicate conclusions regarding 
convection in similar stars based on 1-D model atmo- 
sphere comparisons to optical spectrophotometric data. 
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23.995 


72.69 


81.18 


1.80 


0.10 


1.80 


a L/iVla 


z40zDoo.04o41 


OQ one; 

zo.yyo 


"70 fiO 

/ z.oy 


Q 1 T1 
oi. ( 1 


O TO 

Z. I z 


O OO 

u.uy 


O TO 

Z. 1 z 


a L/iVla 


z4ozDoO.Ooyyo 


OQ onr; 

zo.yyo 


TO fiO 

/ z.oy 


Tn OO 
/ y.zz 


1 RC 

l.Oo 


O OQ 

U.Uo 


1 RC 

l.Oo 


a L/iVia 


z4uoUUz.DoUOZ 


4z.OlcS 


fiO OQ 


t^A KO 

04. Oz 


1 OQ 

l.zo 


O OO 

U.zz 


1 OC 

l.zO 


oU Ljem 


z4ooUUz.Do004 


A Q O/l Q 
4o.y4o 


'70 QO 

/u.oy 


Ki Qn 

01. ©y 


1 1 R 
I.ID 


O OI 

U.zl 


1 1 Q 

l.lo 


oU Ljem 


z4ooUUz.DOiyz 


OU.D4cS 


'70 0*7 
t Z.Zi 


ACt OO 


n K1 
U.Ol 


O 1 Q 

U.lo 


n tiA 
U.04 


oU Ljem 


z4ooUUz.DO 1 o4 


oi.cStty 


'70 CQ 

( z.Oo 


QR QC 
OD.oO 


n '71 
U. ( 1 


O 1 R 
U.lD 


n '7Q 
U. ( o 


Qn (^«™ 
oU Ljem 


z4ozooy.oy 1 oi 


OO. 1 04 


-XAT QQ 

14 (.yy 


QO QO 

oz.yz 


1 QQ 
l.Oo 


U.44 


1 /tn 

1.4U 


Oi OJVia 


z4ozooy.oyooi 


OQ? 

00. zo ( 


1 A'? QO 

14 r .oz 


on Qn 

zy.oy 


o on 

u.yy 


O AC\ 

U.4U 


1 nT 
l.U/ 


a Oivia 


z4ozv5oy .osyu / 


00. / oi 


1 /ifi i^n 

i4D.oy 


Q1 OQ 

oi.zo 


1 11 
1.11 


C\ AO 

U.4z 


1 1 c 

l.lo 


a UiVla 


z40oUUz.oy4z i 


^1 OQO 

/ .zoy 


'7Q QT 
/ O.Oi 


OQ TO 

zo. / y 


o c;o 

u.ou 


O 1 Q 

U. lo 


O c;o 
U.Oz 


oU Liem 


z40v>Uzy.DzU / / 


c;'? oil 


TQ QT 

/ o.o / 


on 1 A 

zy. 14 


O OI 

z.Ul 


O 1 o 
U. lU 


O OO 

z.Uz 


1 C A Ti-in 

lo ivion 


z4oouuz.Dyyou 


c;q 111 
Oo.ili 


VQ /I f; 
( O.40 


OT QC 

z i .yo 


O /I c 
U.4o 


O 1 Q 

U. lo 


O K(\ 

u.ou 


OU *^cin 


z4ooUzy.ozoyo 


c;c 1 Q? 
Oo. iv5 / 


TQ /I c; 
( O.40 


OT TO 

z i . ( y 


1 OO 

i.yz 


O 1 o 
U. lU 


1 OO 

I.yz 


lo ivion 


2452994.73707 


60.383 


73.58 


25.89 


0.45 


0.15 


0.47 


HR 1799 


2453002.72628 


61.657 


73.57 


22.37 


0.73 


0.09 


0.73 


18 Mon 


2453002.73158 


62.172 


73.54 


22.10 


1.05 


0.09 


1.06 


18 Mon 


2452996.75146 


62.477 


73.51 


22.06 


0.25 


0.08 


0.27 


6 Leo 


2452996.75644 


62.863 


73.45 


21.48 


0.26 


0.08 


0.27 


6 Leo 


2452995.76035 


62.946 


73.43 


21.06 


0.33 


0.08 


0.34 


6 Leo 


2452996.81799 


62.952 


71.42 


21.44 


0.33 


0.11 


0.34 


6 Leo 


2452995.81464 


63.326 


71.74 


20.80 


0.27 


0.10 


0.28 


6 Leo 


2452996.81167 


63.339 


71.75 


20.41 


0.27 


0.10 


0.29 


6 Leo 
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TABLE 1 — Continued 



Julia.li Date 


r rOJCCtCQ 


r OSltlOIl 


1/2 
V 








Ojilibration otcir 




Baseline 


Angle 


X 100 


X 100 


X 100 


X 100 






(meters) 


(degrees) 












2452995.80179 


63.845 


72.31 


20.26 


0.27 


0.10 


0.29 


6 Leo 


2452996.78120 


63.942 


72.93 


19.80 


0.18 


0.08 


0.20 


6 Leo 


2453002.76593 


63.957 


72.89 


20.05 


0.59 


0.08 


0.60 


18 Men 


2452996.78597 


63.988 


72.78 


19.85 


0.20 


0.08 


0.22 


6 Leo 


2453002.77184 


63.991 


72.71 


19.88 


0.66 


0.08 


0.66 


18 Men 



TABLE 2 

Mark III 500 nm Visibility Measurements 



Julian Date Projected Position V'^ "'Ko 
Baseline Angle 
(meters) (degrees) 



2447452, 


.956 


2, 


.803 


23, 


.049 


1, 


.024 


0, 


.106 


2447452, 


.970 


2, 


.768 


20, 


.864 


0, 


.900 


0, 


.089 


2447452, 


.986 


2, 


.728 


17, 


.918 


0, 


.944 


0, 


.072 


2447453 


.004 


2, 


.689 


14, 


.543 


0, 


.966 


0, 


.073 


244745.3 


,011 


2, 


.630 


6, 


.448 


1, 


.013 


0, 


.075 


2447453 


.924 


2, 


.877 


27, 


.258 


0, 


.859 


0, 


.058 


2447453 


.936 


2, 


.847 


25, 


.646 


0, 


.911 


0, 


.049 


2447453, 


.949 


2, 


.814 


23, 


.750 


0, 


.984 


0, 


.051 


2447453, 


.962 


2, 


.781 


21. 


.688 


0. 


.904 


0. 


.046 


2447453, 


.981 


2, 


.734 


18. 


.425 


0. 


.940 


0. 


.048 


2447454, 


.921 


5, 


.093 


27. 


.288 


0. 


.894 


0. 


.053 


2447454, 


.957 


4, 


.934 


22. 


.135 


0. 


.780 


0. 


.046 


2447454, 


.969 


4, 


.878 


20. 


.047 


0. 


.804 


0. 


.048 


2447454 


.980 


4, 


.833 


18. 


.154 


0. 


.884 


0, 


.051 


2447454, 


.993 


4, 


.778 


15, 


.552 


0, 


.933 


0, 


.054 


2447455, 


.003 


4, 


.741 


13, 


.490 


0, 


.925 


0, 


.054 


2447455, 


.015 


4, 


.705 


11, 


.116 


0, 


.927 


0, 


.054 


2447455, 


.024 


4, 


.680 


9, 


.063 


0, 


.954 


0, 


.055 


2447455, 


.037 


4, 


.653 


6, 


.184 


0, 


.831 


0, 


.051 


2447455, 


.050 


4, 


.636 


3, 


.183 


0, 


.893 


0, 


.052 


2447455, 


.910 


6, 


.686 


28, 


.266 


0, 


.578 


0, 


.048 


2447455, 


.928 


6, 


.587 


25, 


.979 


0, 


.730 


0, 


.057 


2447455, 


.945 


6, 


.490 


23, 


.593 


0, 


.659 


0, 


.052 


2447455, 


.962 


6, 


.384 


20, 


.683 


0, 


.762 


0, 


.059 


2447455, 


.974 


6, 


.319 


18, 


.661 


0, 


.694 


0, 


.053 


2447455, 


.985 


6, 


.258 


16, 


.546 


0, 


.755 


0, 


.057 


2447455, 


.996 


6, 


.204 


14. 


.363 


0. 


.802 


0. 


.061 


2447456, 


.007 


6, 


.155 


12. 


.079 


0. 


.828 


0. 


.063 


2447456, 


.018 


6 


.114 


9, 


.724 


0. 


.836 


0. 


.064 


2447456, 


.030 


6, 


.079 


7, 


.157 


0. 


.826 


0. 


.063 


2447456, 


.041 


6, 


.055 


4, 


.478 


0. 


.800 


0. 


.061 


2447456, 


.916 


22, 


.193 


27. 


.143 


0. 


.010 


0. 


.002 


2447456, 


.935 


21, 


.836 


24. 


.598 


0. 


.016 


0, 


.003 


2447457, 


.046 


20, 


.215 


2, 


.793 


0, 


.035 


0, 


.009 


2447457, 


.964 


21, 


.220 


19, 


.514 


0, 


.020 


0, 


.002 


2447457, 


.976 


20, 


.984 


17, 


.152 


0, 


.024 


0, 


.003 


2447457, 


.993 


20, 


.703 


13, 


.805 


0, 


.031 


0, 


.002 


2447458 


.008 


20, 


.505 


10, 


.815 


0, 


.030 


0, 


.002 


2447458 


.019 


20, 


.377 


8, 


.311 


0, 


.032 


0, 


.002 


2447458, 


.034 


20, 


.256 


4, 


.843 


0, 


.033 


0, 


.002 


2447458, 


.044 


20, 


.213 


2, 


.678 


0, 


.038 


0, 


.002 


2447458, 


.053 


20, 


.195 


0, 


.435 


0, 


.040 


0, 


.002 


2447460, 


.952 


29, 


.068 


20, 


.149 


0, 


.006 


0, 


.002 


2447460, 


.970 


28, 


.613 


16, 


.858 


0, 


.003 


0, 


.002 



Aufdenberg et al. 



TABLE 3 

Mark III 800 nm Visibility Measurements 



Julian Date Projected Position '^^to 
Baseline Angle 
(meters) (degrees) 



2447452, 


.956 


2, 


.803 


23. 


.048 


1, 


.016 


0. 


.049 


2447452, 


.970 


2, 


.768 


20. 


.863 


0. 


.948 


0, 


.043 


2447452 


.986 


2 


.728 


17, 


.918 


1, 


.009 


0, 


.037 


2447453 


.004 


2 


.689 


14, 


.544 


0, 


.988 


0, 


.035 


2447453 


.022 


2 


.656 


10, 


.718 


0, 


.997 


0, 


.035 


2447453 


.041 


2 


.630 


6 


.447 


1, 


.015 


0, 


.036 


2447453 


.924 


2 


.877 


27, 


.260 


0, 


.964 


0, 


.042 


2447453 


.936 


2 


.847 


25, 


.647 


0, 


.968 


0, 


.029 


2447453 


.949 


2 


.814 


23, 


.750 


1, 


.003 


0, 


.030 


2447453 


.962 


2 


.781 


21, 


.688 


0, 


.963 


0, 


.028 


2447453 


.981 


2 


.734 


18, 


.423 


0, 


.970 


0, 


.027 


2447454, 


.921 


5, 


.093 


27, 


.288 


0, 


.964 


0, 


.027 


2447454, 


.957 


4, 


.934 


22, 


.135 


0, 


.908 


0, 


.026 


2447454, 


.969 


4, 


.878 


20. 


.049 


0, 


.919 


0. 


.026 


2447454, 


.980 


4 


.833 


18. 


.155 


0, 


.967 


0. 


.026 


2447454, 


.993 


4, 


.778 


15. 


.552 


0, 


.986 


0. 


.027 


2447455, 


.003 


4, 


.741 


13. 


.491 


0, 


.978 


0. 


.027 


2447455, 


.015 


4 


.705 


11. 


.115 


0. 


.972 


0. 


.027 


2447455, 


.024 


4 


.680 


9 


.062 


0, 


.977 


0. 


.026 


2447455 


.037 


4 


.653 


6 


.185 


0, 


.933 


0, 


.032 


2447455 


.050 


4, 


.636 


3 


.183 


0, 


.960 


0, 


.027 


244745^ 


.928 


6 


.587 


25, 


.980 


0, 


.881 


0, 


.030 


2447455 


.945 


6 


.490 


23, 


.593 


0, 


.845 


0, 


.030 


2447455 


.962 


6 


.384 


20, 


.683 


0, 


.899 


0, 


.030 


2447455 


.974 


6 


.319 


18, 


.661 


0, 


.856 


0, 


.029 


2447455 


.985 


6 


.258 


16, 


.546 


0, 


.900 


0, 


.030 


2447455, 


.996 


6, 


.204 


14, 


.363 


0, 


.920 


0, 


.030 


2447456 


.007 


6 


.155 


12, 


.078 


0, 


.922 


0, 


.031 


2447456 


.018 


6 


.114 


9 


.724 


0, 


.941 


0, 


.031 


2447456, 


.030 


6, 


.079 


7, 


.157 


0, 


.939 


0, 


.031 


2447456, 


.041 


6, 


.055 


4, 


.478 


0, 


.928 


0, 


.030 


2447438, 


.947 


18, 


.608 


29, 


.260 


0, 


.364 


0, 


.013 


2447438, 


.968 


18, 


.309 


26. 


.817 


0. 


.380 


0, 


.012 


2447438, 


.988 


17, 


.984 


23. 


.980 


0. 


.395 


0. 


.012 


2447439, 


.024 


17, 


.413 


17. 


.976 


0. 


.415 


0, 


.013 


2447439, 


.043 


17, 


.146 


14. 


.289 


0. 


.438 


0. 


.014 


2447456, 


.916 


22, 


.193 


27. 


.143 


0. 


.244 


0, 


.006 


2447456 


.935 


21 


.836 


24. 


.598 


0. 


.250 


0, 


.006 


2447457, 


.046 


20, 


.215 


2 


.793 


0, 


.320 


0, 


.016 


2447457, 


.964 


21 


.220 


19, 


.514 


0, 


.271 


0, 


.006 


2447457, 


.976 


20 


.984 


17, 


.152 


0, 


.290 


0, 


.007 


2447457, 


.993 


20 


.703 


13, 


.805 


0, 


.306 


0, 


.006 


2447458 


.008 


20 


.505 


10, 


.815 


0, 


.312 


0, 


.006 


2447458 


.019 


20, 


.377 


8 


.311 


0, 


.307 


0, 


.007 


2447458 


.034 


20, 


.256 


4, 


.843 


0, 


.335 


0, 


.007 


2447458 


.044 


20 


.213 


2 


.678 


0, 


.332 


0, 


.007 


2447458 


.053 


20 


.195 


0, 


.435 


0, 


.335 


0, 


.007 


2447460 


.952 


29 


.068 


20, 


.149 


0, 


.068 


0, 


.003 


2447460, 


.970 


28, 


.613 


16, 


.858 


0, 


.079 


0, 


.004 


2447478, 


.899 


26, 


.960 


20. 


.833 


0. 


.108 


0, 


.007 


2447478, 


.935 


26, 


.154 


14. 


.128 


0. 


.125 


0. 


.009 


2447478, 


.950 


25, 


.887 


11. 


.006 


0. 


.139 


0, 


.010 


2447478, 


.965 


25, 


.679 


7, 


.698 


0. 


.149 


0. 


.010 


2447484, 


.836 


28, 


.099 


27. 


.661 


0. 


.077 


0. 


.016 


2447484, 


.867 


27, 


.366 


23. 


.477 


0. 


.103 


0, 


.010 


2447484, 


.878 


27, 


.077 


21, 


.631 


0, 


.121 


0, 


.015 


2447484, 


.884 


26 


.926 


20, 


.595 


0, 


.103 


0, 


.011 


2447484, 


.891 


26 


.773 


19, 


.497 


0, 


.117 


0, 


.014 


2447484, 


.897 


26 


.630 


18, 


.398 


0, 


.120 


0, 


.017 


2447484, 


.939 


25 


.806 


9 


.841 


0, 


.143 


0, 


.013 


2447484, 


.951 


25 


.6-59 


7, 


.303 


0, 


.147 


0, 


.015 


2447484, 


.962 


25, 


.552 


4, 


.594 


0, 


.146 


0, 


.014 


2447484, 


.974 


25, 


.493 


1, 


.879 


0, 


.143 


0, 


.012 


2447484, 


.981 


25, 


.482 


0, 


.339 


0, 


.154 


0, 


.013 


2447484, 


.987 


25, 


.486 


-1, 


.146 


0, 


.135 


0, 


.011 


2447484, 


.993 


25, 


.502 


-2, 


.497 


0, 


.124 


0, 


.016 


2447833, 


.958 


24, 


.330 


14 


.949 


0, 


.184 


0, 


.002 


2447833, 


.974 


24, 


.055 


11, 


.678 


0, 


.197 


0, 


.003 


2447833, 


.983 


23, 


.923 


9, 


.706 


0. 


.209 


0, 


.003 


2447833, 


.989 


23, 


.855 


8, 


.511 


0, 


.202 


0, 


.003 


2447833, 


.997 


23 


.767 


6 


.625 


0. 


.206 


0, 


.003 


2447834, 


.002 


23 


.722 


5, 


.442 


0. 


.209 


0, 


.003 


2447834, 


.011 


23, 


.670 


3, 


.558 


0. 


.217 


0, 


.005 
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Julian Date 


Projected 


Position 


1/2 






Baseline 


Angle 






(meters) 


(degrees) 






2447834.024 


23.632 


0.512 


0.214 


0.004 


2447834.029 


23.633 


-0.749 


0.198 


0.004 


2447834.036 


23.648 


-2.353 


0.208 


0.003 


2447834.042 


23.673 


-3.700 


0.203 


0.004 


2447834.049 


23.715 


-5.208 


0.205 


0.003 


2447834.055 


23.766 


-6.598 


0.207 


0.006 


2447834.061 


23.831 


-8.048 


0.182 


0.007 



TABLE 4 

Synthetic Stromgren Indices for Composite SED Components 



Surface Area Weight Tgff(K) log(g) a (6 — j/) ci 



Composite SED A, /3 = 0.6 linear limb darkening 



0.002 


5483 


3.95 


1.25 


0.452 


0.469 


0.027 


5693 


3.95 


1.25 


0.416 


0.446 


0.093 


5893 


3.95 


1.25 


0.384 


0.436 


0.125 


6074 


3.95 


1.25 


0.356 


0.445 


0.120 


6241 


3.95 


1.25 


0.332 


0.464 


0.108 


6395 


3.95 


1.25 


0.311 


0.491 


0.114 


6538 


3.95 


1.25 


0.291 


0.521 


0.134 


6672 


3.95 


1.25 


0.274 


0.554 


0.139 


6799 


3.95 


1.25 


0.257 


0.586 


0.106 


6920 


3.95 


1.25 


0.242 


0.620 


0.029 


7034 


3.95 


1.25 


0.227 


0.652 


0.002 


7145 


3.95 


1.25 


0.213 


0.684 



Composite SED B, /3 = 0.5 linear limb darkening 



0.002 


5539 


3.95 


1.25 


0.442 


0.461 


0.027 


5751 


3.95 


1.25 


0.406 


0.444 


0.093 


5953 


3.95 


1.25 


0.375 


0.437 


0.125 


6136 


3.95 


1.25 


0.347 


0.451 


0.120 


6304 


3.95 


1.25 


0.323 


0.474 


0.108 


6460 


3.95 


1.25 


0.302 


0.504 


0.114 


6604 


3.95 


1.25 


0.283 


0.536 


0.134 


6740 


3.95 


1.25 


0.265 


0.571 


0.139 


6869 


3.95 


1.25 


0.248 


0.606 


0.106 


6990 


3.95 


1.25 


0.233 


0.639 


0.029 


7106 


3.95 


1.25 


0.218 


0.672 


0.002 


7217 


3.95 


1.25 


0.205 


0.704 



TABLE 5 
Fit Angular Diameters 



Model 


Parameters 


Mark III 


Mark III 


VLTI/ VINCI 






500 nm 


800 nm 


2.2^jm 






(mas) 


(mas) 


(mas) 



"A" PHOENIX 
"B" PHOENIX 
"CI" CO"'BOLD 
"C2" CO^BOLD 
"D" ATLAS 12 
"E" ATLAS 12 
"F" ATLAS 12 

"U" Uniform Disk 



Toff =6530 K, log(g) = 3.95, a = 1.25 
Teff =6530 K, \og(g) = 3.95, a = 0.5 
PHOENIX Tcff =6500 K, \og{g) = 4.0, global mean structure 
PHOENIX Tcff =6500 K, \og{g) = 4.0, 12 weighted structures 

Te£f=6530 K, \og{g) = 3.95, a = 1.25, no overshooting 
Tefi=6530 K, \og{g) = 3.95, a = 1.25, 100% overshooting 
refi=6530 K, log(g) = 3.95, a = 1.25, 50% overshooting 

y2 = |2Ji(7r0uDS/A)/(7reuDB/A)|^ 



5.449±0.012 
5.468±0.012 
5.387±0.011 
5.409±0.011 
5.426±0.011 
5.368±0.011 
5.400±0.011 



5.416±0.005 
5.415±0.005 
5.388±0.005 
5.393±0.005 
5.404±0.005 
5.383±0.005 
5.398±0.005 



5.411±0.006 
5.410±0.006 
5.402±0.006 
5.405±0.006 
5.404it0.006 
5.394±0.006 
5.402±0.006 



5.044±0.010 5.208±0.005 5.302±0.005 
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TABLE 6 
Synthetic Stromgren Indices 



Model 


Tcff (K) 






Additional Parameters 


(o - y) 


ci 


ATLAS 12 "E" 


6530 


3.95 


a 


= 1.25, overshooting 


0.321 


0.378 


ATLAS 12 "F" 


6530 


3.95 


a 


= 1.25, 50% overshooting 


0.311 


0.402 


ATLAS 12 "D" 


6530 


3.95 


a 


= 1.25, no overshooting 


0.302 


0.437 


PHOENIX "A" 


6530 


3.95 


a 


= 1.25 


0.292 


0.519 


PHOENIX "B" 


6530 


3.95 


a 


= 0.50 


0.282 


0.547 


Composite SED A 


6466 


3.95 


/3 


= 0.60 linear limb darkening 


0.293 


0.539 


Composite SED B 


6531 


3.95 


/3 


= 0.50 linear limb darkening 


0.283 


0.553 



" I nterstellar reddening towards Procyon is assumed to be negligibl e: (b — y) = (b — y)^, 
iCirawford fc Barnes 1970: Nordstrom ct al. 200ll and ci = 0.532 (Cr awford fc BarneiilT97B) . 



cq. Observed values: {b — y) = 0.272 



TABLE 7 

Fundamental Parameters for Procyon A 



Parameter 



Symbol 



Value 



Reference 



Limb-darkened angular diameter (mas) 
Bolometric flux (erg cm~^ ^~^) 
Effective temperature (K) 
Parallax (mas) 
Radius (Rq) 
Mass (Mq) 

Surface gravity (cm s~-^) 



^LD 5.404±0.031 

(17.86 ± 0.89) X lO-'' 

Tcff 6543±84 

TT 285.93±0.88 

R 2.031±0.013 

M 1.42±0.04 

log(g) 3.975±0.013 



Table O §E1 

9ld and 
Pe^^mar^^al 

Girard et al. (2000 
R and M 
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X [Mm] 

Fig. 1. — Typical map of the emergent white hght intensity in vertical direction encountered during the temporal evolution of the grey 
CQ^BOLD model. The black line depicts an isophote at 3.95 X 10^'' erg cm~^ sr~^ s^^ approximately corresponding to the image's median 
intensity (see Figure 1^. 
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Fig. 2. — Probability density function (PDF), the probability per differential intensity interval with the units of inverse intensity, of the 
emergent white light intensity of the grey CO^BOLD model. The PDF was derived from a sequence of intensity maps like the one shown in 
Figure^ Twelve intensity intervals are demarcated by the dashed lines. The surface area weights for these intervals are given in Table Bl 
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Fig. 3. — (a) The squared visibility data from Tables (red), |^ (purple) and |21 (green) as a function of projected baseline compared 
to the synthetic visibilities (from eq.(2)) from the 3-D CO^BOLD + PHOENIX model "C2" with an angular diameter of 5.403 mas. (b) The 
deviations of each data set from the model, (c) Histograms of these deviations for each data set and for all the data combined. 
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Fig. 4. — A comparison of the best fit angular diameters ^ld 500 nm, 800 nm (Mark HI data), and 2.2 ^m (VLTI data) from Table|Sl 
for the seven atmosphere models "A" to "F" and one uniform disk model "U" . The uniform disk angular diameters indicate that Procyon 
is more limb darkened at shorter wavelengths as expected. The atmosphere model fits are more dispersed at shorter wavelengths. The 
fits indicate that f-D models without overshooting are too limb darkened at 500 nm, while f-D model "E" with f00% overshooting is not 
limb darkened enough. Model "F" with 50% overshooting yield the same angular diameter (within la) at all three wavelengths. The 3-D 
CQ^BOLD + PHOENIX models "Cf" and "C2" yield angular diameters slightly less consistent than model "F", but with no free parameters 
for convection. The dashed line indicates the angular diameter of the model in Figure 1^ 
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Fig. 5. — Comparisons between synthetic spectral energy distributions and spectrophotometric measurements at ultraviolet and visual 
wavelengths. The models and data are binned to 2 nm resolution in the UV for clarity. The models are binned to 5 nm in the visual to 
match the resolution of the spectrophotometry. The synthetic SEDs are scaled in absolute flux using an angular diameter of 5.40 mas, 
a value consistent with our best interferometric estimate (see Table 0. (a) Composite SED model B (green), which incorporates a 3-D 
model for the distribution of surface intensities due to granulation, better represents the UV flux distribution than a single T^g component 
1-D model (blue, ATLAS 12 model "E"), particularly at wavelengths below 160 nm, where all five single T^ff component, 1-D atmosphere 
models fail, (b) Comp osite SED model B (green) yields Stromgren v- and b-band fluxes consistent with the observed spectrophotometry 
IGlushneva et alJ1992i) and larger than the 1-T(,fi component model (blue, ATLAS 12 model "E"). These two models are in closer agreement 
in the u- and jy-bands. 
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Fig. 6. — A comparison of the synthetic (b — y) and ci Stromgren indices from TablelSlrelative to the observed indices ICrawford fc Earned 
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Fig. 7. — A comparison of six model atmosphere temperature on a mean Rosseland optical depth scale: PHOENIX models "A" and "B" 
(dash-dot), ATLAS 12 models "D", "E", "F" (dot), and CO^BOLD model "CI" (solid). The degree of limb darkening, particularly at 500 nm 
where the opacity is relatively low compared to 800 nm and where the intensity contrast as a function of temperature is relatively high 
compared to 2.2 fim, is sensitive to different convection treatments which affect the temperature gradient in the interval tRoss = 1 to ~ 3. 
The model with the steepest gradient, PHOENIX "B" with a = 0.5 and no overshooting, has the strongest limb darkening. ATLAS 12 model 
"E" with 100% overshooting has the shallowest gradient and the weakest limb darkening. ATLAS 12 model "F" and CO BOLD model "CI" 
have similar gradients and provide an intermediate degree of limb darkening which is consistent with the interferometric observations. 
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0.01 0.10 1.00 10.00 100.00 

Rosselond Mean Opticol Depth 

Fig. 8. — A comparison of the six model atmosphere radiative flux structures on a mean Rosseland optical depth scale: PHOENIX models 
"A" and "B" (dash-dot), ATLAS 12 models "D", "E", "F" (dot), and CD^BDLD model "CI" (solid). The 1-D overshooting models "E" and 
"F" and the temporally and spatially averaged 3-D model "CI" all have a significant fraction (5% to 15%) of their total flux in convection 
at an optical depth of unity, while the other 1-D models are fully radiative at this depth. 



